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ABSTRACT 

O , We investigate the flux density variability and changes in the parsec-scale radio struc- 

ture of the flat spectrum radio quasar PKS 1510-089. This source was target of multi- 
epoch Very Long Baseline Interferometer (VLBI) and Space- VLBI observations at 4.8, 
8.4 and 22 GHz carried out between 1999 and 2001. The comparison of the parsec- 
scale structure observed at different epochs shows the presence of a non-stationary 
jet feature moving with a superluminal apparent velocity of 16.2c ±0.7c. Over three 
epochs at 8.4 GHz during this period the core flux density varies of about 50%, while 
the scatter in the jet flux density is within 10%. The polarization percentage of both 
core and jet components significantly change from 2 to 9 per cent, while the polariza- 
tion angle of the core shows an abrupt change of about 90 degrees becoming roughly 
perpendicular to the jet direction, consistent with a change in the opacity. To complete 
the picture of the physical processes at work, we complemented our observations with 
multi-epoch Very Long Baseline Array (VLBA) data at 15 GHz from the MOJAVE 
programme spanning a time baseline from 1995 to 2010. Since 1995 jet components 
are ejected roughly once per year with the same position angle and an apparent speed 
between 15c and 20c, indicating that no jet precession is taking place on a timescale 
longer than a decade in our frame. The variability of the total intensity flux density 
together with variations in the polarization properties may be explained assuming ei- 
ther a change between the optically-thick and -thin regimes produced by a shock that 
varies the opacity, or a highly ordered magnetic field produced by the compression of 
the relativistic plasma by a shock propagating along the jet. Taking into account the 
high 7-ray emission observed from this source by the AGILE and Fermi satellites we 
investigated the connection between the radio and 7-ray activity during 2007-2010. 
Multi-wavelength flux and polarization observations suggest that during some 7-ray 
flaring episodes the emission at high (7-ray) and low (radio) energies has origin in the 
same region. 

Key words: polarization - galaxies quasars: individual (PKS 1510-089) - radio con- 
tinuum: general - radiation mechanisms: non-thermal 
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1 INTRODUCTION 

Powerful radio emission is rare in the population of active 
galactic nuclei (AGN). This phenomenon, only found in 
about 10% of the AGNs, is associated with the presence of 
relativistic particles produced in the central regions of the 
AGN, and then channelled through the jets out to distances 

* E-mail: orienti@ira.inaf.it 



of hundreds of kpc up to a few Mpc. Such relativistic 
particles are responsible for the synchrotron radiation at 
the origin of the radio emission. Furthermore, they may 
scatter low energy photons to high energy bands by inverse 
Compton (IC) processes, suggesting a connection between 
radio and high energy 7-ray emission. Indeed, all the 7-ray 
AGN detected by EGRET are radio-loud sources, strongly 
supporting this scenario. No radio-quiet AGNs were firmly 
detected in 7-rays by the Large Area Telescope (LAT) on 
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board the Fermi satellite during the first year of operation, 
although possible association for the LAT sources with 
radio-quiet A GNs are tentatively proposed in the first LAT 
AGN Catalog l|Abdo et alj|2010ah . As Fermi-L AT continues 
to collect data and its sensitivity increases, some of these 
associations could be confirmed. 

The AGN dominating the 7-ray sky is the blazar popula- 
tion, comprising flat spectrum radio quasars (FSRQ) and 
BL Lac objects. These sources are characterized by the 
presence of a compact radio core, apparent superluminal 
jet speed, extreme flux density variability in all bands, 
and high fraction of polarized optical and radio emission. 
Their observational properties are interpreted as the result 
of severe beaming effects due to the orientation of the 
relativistic jet at very small angles to the line of sight. 
Thanks to the episodes of enhanced luminosity across 
the entire electromagnetic spectrum, it is possible to set 
constraints on the physical properties of the region along 
the jet responsible for the emission at the various wave- 
lengths. Radio monitoring of EGRET sources suggested 
that the highest levels of 7-ray emission is observed close i n 
time to radio flares (e.g. lLahteenmaki fc Valtaoial 12003') . 
and connected with the emission of a new jet component 
ijjorstad et alJl200ll ). Both pieces of evidence give support 
to the idea that the strongest 7-ray emission is strictly 
related to a sho ck in the jet that produ ces the synchrotron 
radio flare (e.g. iMarscher fc Geanll985l ). It must be noted 
that the EGRET data were sparse, with large uncertainties 
and selection effects, precluding a clear temporal correlation 
between 7-ray activity and the emission at lower frequen- 
cies. The advent of the AGILE and Fermi 7-ray satellites 
allowed us to test the results obtained in the EGRET era, 
providing details on the connection between 7-ray and radio 
emission to find out the mechanisms responsible for the 
high-energy emission. Correlation studies of the first three 
months of Fermi-ljAT data and quasi-simultaneous Very 
Long Baseline Interferometry (VLBI) observations showed 
that the ga mma-ray emitting b lazars have faster apparent 
jet spe eds dLister et al.l l2009al ). wider apparent opening 
angles (|Pushkarev et al1l2009h. and hig her variability and 
Doppler factor ( Savolainen et all |2010| ). with respect to 



blazars with weak 7-ray emission. Variability studies from 
radio to high energy bands give important insight to locate 
the site of flares and infer the physical conditions of the 
jet. In particular, using the most recent 7-ray as well as 
radio polarimetric data, a likely close connection between 
high 7-ray states and the activity i n the pc-sca l e cor e 
region has been investigated (see e.g. lAgudo et al.ll201lT ). 
In addition, Pushkarev et al. (2010) found for a sample of 
186 sources a non-zero time delay between radio emission 
measured by pc-scale observations at 15 GHz and 7-ray 
radiation detected by Fermi-LAT, suggesting that the delay 
is most likely connected with synchrotron opacity in the 
core region. 

Among blazars, PKS 1510-089 is an ideal target to inves- 
tigate the location of the emitting region and the physical 
processes occurring in relativistic jets. This object is a 
flat s pectrum radio quasar at z = 0.361 (| Thompson et all 
1990) with highly polarized optical emission. The syn- 
chrotron emission has its peak in the IR band, while the 
inverse Compton component peaks in the 7-ray regime. 
High flux density variability is observed throughout the 



electromagnetic spectrum, from the radio to the high 
energy bands. In particular, episodes of hig h 7-ray activity 
were r ecently detected by both Fermi-LAT (ICutini fc Havsl 
120091 ; ICiprini fc Corbell 120091 ; iTramacerd |200Sh and the 
Gamma Ray Imagi n g Detector (GRID) on board AGILE 
dStriani et all l20ld; IVercellone et ail 120091 ; iPucella et ail 
120091 ; iD'Ammando et alj l2009al . l200Sh . Description of 
these events, together with their possible co nnection with 
multiwavelengths emission can be found in IPucella et al] 
J200S|). ID'Ammando etall (l2009bl). IMarscher et al.l i|2010t ). 
lAbdo et all i2010bj ). and lD'Ammando et al.l (|201ll ). 
In the radio band the emission is dominated by the 
core component. Multi-epoch parsec-scale observations 
revealed highly su perluminal knots wi t h apparent velocity 
exceeding 2 c (e.g. iHoman et aill200ll; Ijorstad et al.l 120051 : 



iLister et all l2009bl ; IMarscher et al.l |2010h . ejected along 
the north-west direction at an angle of about -30° with 
respect to the core. On arcsecond scale the jet structure 
is oriented in the opposite direction, indicating a severe 
misalignment of almost 180° between the pc- and kpc-scale 
jet. Misalignment between pc- and kpc-scale structure 
has been frequently observed i n core-dominated sources 
(e.g. IPearson fc Readheadl Il988l ). However, misalignment 
larger than 110° a s those found in 0954+556 and 1652+398 
l|Lister et al.ll200~lh is very rare. The extraordinary case o f 
PKS 1510-089 has been described bv lHoman et all (|2002aT ) 
assuming a small change of about 12°-24° in the intrinsic jet 
direction, that appears amplified due to projection effects, 
providing a simple explanation for the observed morphol- 
ogy. Both the misalignment and the highly superluminal 
jet speed indicate that the jet axis of PKS 1510-089 forms a 
very small angle of a few degrees to the line of sight. Such an 
extreme orientation enhances beaming effects making this 
source a good target to investigate the possible connection 
between 7-ray flares and the ejection of jet components. 
In this paper we present new results of proprietary multi- 
epoch polarimetric VLBI and Space- VLBI observations 
of PKS 1510-089 carried out at 4.8, 8.4 and 22 GHz 
between 1999 and 2001 and not yet published, with the 
aim of studying changes in the source structure and 
investigating their possible connection with other physical 
properties, such as flux density variability, spectral index 
distribution and polarization properties. We then compare 
our results with multi-epoch Very long Baseline Array 
(VLBA) data at 15 GHz from the MOJAVE (Monitoring 
Of Jets in Ac tive galactic nuclei w ith VLBA Experiments) 
programing]] (|Lister et al. I [2009c) spanning a larger time 
interval, between 1995 and 2010. The addition of high 
energy information from AGILE and Fermi data gives us 
a clue to connect 7-ray emission and radio properties, like 
flux density variability and changes in the pc-scale radio 
structure for the period 2007-2010. 

Throughout this paper, we assume the following cos- 
mology: Ho = 71 km s" 1 Mpc~\ Q M = 0.27 and Ov = 0.73, 
in a flat Universe. At the redshift of the target 1 arcsec = 



1 All the 15-GHz data presented in this paper are from the MO- 
JAVE programme. The MOJAVE data archive is maintained at 
|http : / /www, physics . purdue. edu / M O J AVE 
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Figure 1. Plot of amplitude vs. projected baseline length at 4.8 
GHz (top), and the uv coverage (bottom) of Space- VLBI (VSOP) 
observations. 



5.007 kpc. The spectral index a is defined as S(y) oc v a . 
2 RADIO DATA 

2.1 VLBI and Space- VLBI observations 

PKS 1510-089 was target of Space- VLBI observations. It was 
observed at 4.8 GHz (C band) by VLBA+HALCA on 1999 
August 11 and on 2000 May 13 in single polarization mode. 
In each observing run the target was observed for about 6 
hours. Examples of the uv coverage and amplitude vs base- 
line length are presented in Fig. [1] 

The target source was observed with the VLBA with the 
addition of the Effelsberg telescope at 8.4 GHz (X band) 
and 22.2 GHz (K band), in full polarization mode with 
a bandwidth of 32 MHz at 128 Mbps. To study possible 
changes in the source structure, subsequent VLBA observa- 
tions at 8.4 GHz were carried out in full polarization with a 
recording bandwidth of 32 MHz at 128 Mbps. The correla- 
tion was performed at the VLBA correlator in Socorro and 
the data reduction was carried out with the NRAO AIPS 



package. In the X band the instrumental polarization was 
removed by using the AIPS task PCAL. The absolute orien- 
tation of the electric vector of the calibrator B1749+096 was 
then compared with the VLA/VLBA polarization calibra- 
tion database to derive the corrections. The values derived 
are in good agreement 5°). 

Final images for each epoch and at each frequency 
were produced after a number of phase self-calibration 
iterations. Amplitude self-calibration was applied at the 
end of the process using a solution interval longer than 
the scan length, to remove residual systematic errors and 
to fine tune the flux density scale, but not to force the 
individual data points to follow the model. The resolution 
at the various frequencies are almost comparable (Table [1} 
due to the similar wti-coverage reached with the different 
interferometers used. At 8.4 GHz, besides the total inten- 
sity (I), images in the Stokes' U and Q parameters were 
produced with the final fully calibrated datasets. Final 
VLBI images at 4.8, 8.4 and 22.2 GHz are shown in Fig. [2] 
For the observations carried out on 1999 January 11, we 
produced also a low-resolution image at 22.2 GHz using 
the same uv range, restoring beam and image sampling of 
the 8.4 GHz data in order to produce the spectral index 
image, which is presented in Fig. [3] superimposed on the 
22.2 GHz contours obtained from the low-resolution image. 
Information on the VLBI observations is reported in TableQ] 



2.2 MO J AVE data 

To study the radio variability of PKS 1510-089, we com- 
pared our observations with multi-epoch 15-GHz (U band) 
VLBA data from the MOJAVE programme spanning a 
time interval from 1995 July 28 to 2010 December 23. The 
typical resolution is about 1.4x0.5 mas. For each of the 51 
epochs analy s ed, we imported the calibrated tt-y-datasets 
l|Lister et al. I |2009ch into the NRAO AIPS package and 
performed a few phase-only self-calibration iterations before 
producing the final total intensity images w hich resulted 
to be f ully consistent with those reported by iLister et al. I 
|2009d\ For those datasets in full polarization mode, we 
produced also Stokes' U and Q images. The rms noise level 
measured on the image plane is in the range of 0.15 and 0.3 
mjy/beam. Total intensity images concerning the observing 
epochs betw e en July 1995 and July 2008 are published in 
ILister et al. I l|2009ci ). An image of the source structure in 
September 2010 is presented in Fig. [4] as an example. 



2.3 Data analysis 

The flux density and deconvolved angular size of each 
source component were measured by means of the AIPS 
task JMFIT, which performs a Gaussian fit to the source 
components on the image plane. For extended components 
like the low-surface brightness diffuse jet visible at 8.4 GHz, 
the flux density was derived by means of TVSTAT, which 
performs an aperture integration over a selected region on 
the image plane. When the polarization information was 
available, we derived the polarization parameters, such as 
the linearly polarized flux density, the fractional polariza- 
tion, and the (integrated) position angle of the electric 
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Figure 2. VLBI images of the source PKS 1510-089. Upper panel: total intensity VSOP images at 5 GHz and VLBI image at 22 GHz. 
Middle panel: total intensity VLBA images at 8.4 GHz. Bottom panel: fractional polarization images at 8.4 GHz. On each image we 
provide the telescope array, the observing band and the observing date; the peak flux density in mjy/beam and the first contour (f.c.) 
intensity in mJy/bcam, which corresponds to three times the off-source noise level. Contour levels increase by a factor 2 in the upper 
and middle panels, and a factor 4 in the bottom panel. The restoring beam is plotted on the bottom left corner. For the datascts with 
polarization information, the vectors superimposed on the I contours show the position angle of the E vector, where one-millimeter 
length corresponds to 2.5 mjy/beam. The grcyscale represents the fractional polarization. 
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Figure 3. Spectral index distribution between 8.4 and 22.2 GHz 
across the target source PKS 1510-089 superimposed on the 22.2 
GHz contours convolved to the 8.4-GHz beam. On the image we 
provide the peak flux density in mjy/beam and the first con- 
tour (fx.) intensity in mjy/beam corresponding to 3 times the 
off-source noise (lcr). Contour levels increase by a factor 2. The 
restoring beam is plotted on the bottom left corner. 
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Figure 5. VLA A-array image at 8.4 GHz of PKS 1510-089 pro- 
duced from archival data (project AF376) and calibrated using 
the standard procedure developed in the NRAO AIPS package. 
On the image we provide the observing band, the peak flux 
density in mjy/beam and the first contour (f.c.) intensity in 
mjy/beam, which corresponds to three times the off-source noise 
level. Contour levels increase by a factor 2. The restoring beam 
is plotted on the bottom left corner. 



PKS 1510-089 VLBA U band 29/09/2010 



£ < - 



peak= 1975.8; f.c.= 0.7 (mjy/beam) 




MilliARC SEC 



Figure 4. An example of the VLBA image at 15 GHz relative 
to 2010 September 29. On the image we provide the observing 
band, the peak flux density in mjy/beam and the first contour 
(f.c.) intensity in mjy/beam, which corresponds to three times 
the off-source noise level. Contour levels increase by a factor 2. 
The restoring beam is plotted on the bottom left corner. 



vector E. In Table [2] we report the source parameters 
concerning our proprietary observations at 4.8, 8.4, 22 GHz 
and listed in Table Q] 

In addition to the data analysis performed on the image 
plane, a further modelfitting with Gaussian components 
was carried out to the visibility data at each epoch using the 
modelfitting option in Difmap. This approach is preferable 



in the case we want to derive small variations in the source 
structure. Indeed, the analysis of the visibility data allows 
us to use the full resolution capability of the interferometer, 
without the dependence of the beam and sampling used 
in producing the image. Furthermore, it provides more 
accurate fit of unresolved structures close to the core 
component, as the case of new jet features, thus yielding 
to a better determination of the position of each component. 



3 RESULTS 
3.1 Morphology 

The radio emission of PKS 1510-089 is dominated by the 
bright core component from which a pc-scale jet emerges 
at an angle of about -30° (Fig. [2}. The pc-scale jet is not 
straight, but it bends at about 2 mas (10 pc) from the 
core, and its emission progressively fades until 5 mas (25 
pc) where it goes below the detectio n limit, in agreement 
with images at oth e r frequencies (e.g. Homan et alj [2002a ; 
I Jorstad et all 120051 ; iLister et al~1 l2009cT l. Furthermore, an 
additional component almost perpendicular to the jet axis, 
visible in several observing epochs and labelled P in Fig. 
[5] is present at about 2 mas West from the core. This 
component is present only in the 15-GHz images with the 
highest signal-to-noise level precluding a reliable discussion 
on its nature. 

In our deep VLBA observations at 8.4 GHz (Fig. [2} a 
low-surface brightness diffuse jet is visible up to about 25 
mas (125 pc) from the core. 

From the spectral index image (Fig. |3j produced comparing 
the simultaneous observations at 8.4 and 22.2 GHz we find 



6 M. Orienti et al. 

Table 1. VLBI observations of PKS 1510-089. 



Freq 


Code 


Obs. date 


Dur. 


rms 


Beam 


Array 




GHz 






h 


mjy/b 


mas X mas 






4.8 


W004 


1999 Aug 11 


6 


0.97 


1.80x0.90 


VLBAH 


-HALCA 


4.8 


W004 


2000 May 13 


6 


0.57 


1.80x0.90 


VLBAh 


-HALCA 


8.4 


BV027 


1999 Jan 11 


1 


0.90 


2.29x0.88 


VLBA-i 


-EF 


8.4 


BV042 


2000 Dec 7 


1.5 


0.23 


2.44x1.00 


VLBA 




8.4 


BV042 


2001 May 17 


1.5 


0.27 


2.36x0.99 


VLBA 




22.2 


BV027 


1999 Jan 11 


1.5 


0.47 


1.70x0.49 


VLBAH 


-EF 



Table 2. VLBI observational parameters of PKS 1510-089 for the observing epochs reported in Table 1. 



Comp I/ obs 5 cp i 5'cp2 S cp 3 Pol cp i Pol cp 2 Pol ep 3 Xcpl Xop2 Xcp3 

GHz mjy mjy mjy mjy (%) mjy (%) mjy (%) deg deg deg 



C 4.8 1504 1161 - - - - 

C 8.4 1120 807 1560 98 (8.6%) 23 (2.8%) 26 (1.6%) 140 110 16 

C 22.2 1139 - - - 



J 


4.8 


626 


644 










J 


8.4 


236 


261 


244 15 (6.3%) 13 (5.0%) 7 (2.9%) 


70 


77 


100 


.1 


22.2 


143 













Diffuse J 8.4 - 34 70 



that the core component has a slightly inverted spectrum 
(q = -0.02), while in the jet the spectral index is 0.5 close 
to the core and steepens moving outwards. 
This initial part of the jet (0.3", i.e. 1.5 kpc) is also visible 
in the VLA image at 8.4 GHz (Fig. as the "northern 
component" CI, in the opposite direction of the kpc-scale 
structure (labelled S). The large scale structure is roughly 
collimated up to a distance of about 2 arcsecond (10 kpc) 
from the core. Then it slightly bends to the West, forming 
an extended low-surface brightness lobe-like feature. Two 
compact regions (SI and S2 in Fig. 0, likely jet knots, are 
located at 1.0" (~ 5 kpc) and 2.7" (~13 kpc) from the core. 

3.2 Flux density and polarization 

The analysis of the radio lightcurves of PKS 1510-089 
shows strong flux den sity variability, wh e re low states 
alternate with flares (jVenturi et al.l l200ll ; Ijorstad et al.l 
l200ll ; iHughes et alj fl992l ) . Sometimes additional episodes 
of abrupt flux density increase have been registered. These 
studies are based on single-dish observations where it is not 
possible to disentangle the contribution of the jet from the 
one arising from the core, and changes in the parsec-scale 
structure may be washed out by the contribution of the 
stationary components. The availability of high-resolution 
VLBI data allows us to separate the core and jet emission 
and thus to determine how the flux density changes in 
both components. Comparing our three-epoch data at 8.4 
GHz (Table [1} we find that the core flux density varies of 



about 50% between 1999 and 2001, while the scatter in 
the jet flux density is within 10%. We compare the flux 
density lightcurve at 8.4 GHz with the 15-GHz VLBA 
data from the MOJAVE programme l|Lister et al. I |2009(J ) 
obtained between 1998 and 2001. In Figs. [B^,b we report 
the flux density of the core and jet components. The core 
flux density at 15 GHz clearly shows strong variability 
reaching a maximum in October 1999, when it has almost 
doubled its value, then followed by a clear decrease lasting 
till the end of 2001 (Fig. After this period the core 
flux density increases again, in agreement with the trend 
found at 8.4 GHz, but the observations are quite sparse 
precluding a more detailed study of the core variability. 
On the other hand, the jet flux density at 15 GHz varies 
about 25% till the end of 2000, while the scatter increases 
in 2001 (Fig. [Sp). The jet flux density has more moderate 
variability than the core. Despite the lack of simultaneous 
multifrequency observations, we can confirm a steep overall 
spectral index for the jet (a —0.5 - 0.7) between 8.4 and 
15 GHz. On the other hand a characteristic value of the 
spectral index of the core could not be determined due to 
its strong and irregular variability associated with changes 
in the opacity. A realistic value for the core spectral index 
could be derived only for the observations carried out on 
January 1999 when simultaneous observations at 8.4 and 
22 GHz are available, providing a flat spectral index (Fig. 

E}. 

The polarization properties of both core (Fig. [BJ:,e) and 
jet (Fig. [6ji,f) components are quite different. In the three 
observing epochs at 8.4 GHz we find that both the core 
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Figure 6. Total intensity and polarization properties at 8.4 and 
15 GHz of PKS 1510-089. Core parameters are marked by a dia- 
mond at 15 GHz, and by a cross at 8.4 GHz. Jet parameters are 
marked by a plus sign at 15 GHz and by an asterisk at 8.4 GHz. 
Panel a: total intensity flux density of the core; panel b total 
intensity flux density of the jet; panel c: fractional polarization 
of the core component; panel d: fractional polarization of the jet 
component; panel e: position angle of the electric vector (EVPA) 
of the core component; panel f: EVPA of the jet component. Un- 
certainties are within the symbols. 



and the jet are polarized with a fractional polarization 
varying between 1.6 and 8.6 per cent in the core region, 
and about 2.9 and 6.3 per cent in the jet (Table [21 bottom 
panel of Fig. [2] and Fig. [6t,d). At 15 GHz the polarization 
percentage of the core is very variable, between 1.6 and 
5.5 per cent (Fig. [BJ;), and shows no evident correlation 
with the total intensity flux density. On the other hand, 
the polarization percentage of the jet at 15 GHz (Fig. 
[6ji) is roughly constant around 4.5 per cent. In the core 
component the position angle of the electric vector E 
(EVPA) at 8.4 GHz is 140° and 110° during the first two 
epochs, then it changes abruptly in May 2001, when it is 
16°. A similar behaviour is found at 15 GHz, where the 
EVPA is between 168 and 130 degrees from 1999 to 2000, 
i.e. roughly parallel to the jet axis, while in June 2001 it 
is 32 degrees (Fig. [6^), becoming perpendicular to the jet 
direction. On the other hand, the jet EVPA does not show 
such a large variation, being between 80 and 100 degrees, 
with the exception of September 2000 when it turned 
out to be 120 degrees. We note that only the integrated 
rotation measure (RM) is available for this source (-15±1 
rad m~ 2 , Simard-Normandin et al. 1984). The corrections 
to the observed EVPA to obtain the intrinsic orientation 
are therefore negligible. 



3.3 Proper motion 

The multi-epoch analysis of the pc-scale morphology of 
PKS 1510-089 shows a considerable evolution of the source 
structure: jet components emerge from the core at different 
times and their changes can be followed by comparing 
observations carried out after short time intervals. With the 
aim of characterizing variations in the source structure we 
modelfitted the visibility data at each epoch (see Section 
2.3). Direct comparison of models obtained independently 
at each epoch is not the b est approach to detect small 
changes (| Conway et alj[l993 ). For this reason, we produced 
a zero-order model consisting of 3 elliptical Gaussian com- 
ponents, which was used as the initial model in modelfitting 
the visibility data of each observing epoch. Errors Ar 
associated with the component position were estimated by 
means of: 



Ar = a/(Sp/rms) 



(1) 



where a is the component deconvolved major axis, 5 P is its 
peak flux densit y and rms is the la noise le vel measured on 
the image plane jPolatidis &: Conw ay 2003). In the case the 
errors estimated by Eq. [1] are unreliably small, we assume 
a more conservative value for Ar that is 10% of the beam. 
The data points are then fitted by a linear model that 
minimizes the chi-square error statistic. 
The linear fit on the three epochs of 8.4-GHz data (Table 
[T] and Fig. [5J) provides an angular separation rate of 
0.990±0.040 mas/yr that corresponds to /3j = 16.2 ± 0.7. 
From the linear back extrapolation fit, and considering the 
uncertainties on the fit parameters, we estimate the time 
of zero separation To between the jet component and the 
core, that results to be T ,j =1997.42±0.12. The accuracy 
of the fit has been tested also by comparing the component 
separation derived by modelfitting the visibility of the two 
epochs 4.8-GHz Space- VLBI data and the three epochs 
at 15 GHz available from the MOJAVE project in the 
same time interval (left panel in Fig. [7]), all with a similar 
resolution. 

To extend the analysis of the knot separation speed till the 
end of 2010, we modelfitted the data at 15 GHz from the 
MOJAVE programme. Results from the analysis of data 
between 1 995 and 2008 were a lready published in previou s 
works by iHoman et ail (|200lh and iLister et ail (|2009bh . 
where the separation velocities found for the detected knots 
range between 15c and 20(0- Since 2007, when also Fermi 
and AGILE satellites could provide crucial information to 
relate the radio and high-energy 7-ray emission, we could 
follow the evolution of three additional knots, labelled Nl, 
N2 and N3 in Fig. [4] For Nl, a l a separation speed of 675 
jtias/yr (i.e. 11.6c) was reported bv lAbdo et al.1 (|2010bT ) using 
data between September 2008 and December 2009. In our 
analysis we expand the time baseline to December 2010. We 
determine the angular separation speeds of these three new 
knots from the core, which we assume stationary: we derive 
1.060±0.056 mas/yr, 1.102±0.113 mas/yr, and 1.041±0.250 
mas/yr for Nl, N2 and N3 respectively, which correspond 
to /9 N i = 17.3 ± 0.9, /3n2 = 18.0 ± 1.9, /3 N3 = 17.0 ± 4.0. 

2 In the Appendix, we re-modelfitted these datasets in order to 
compare results obtained with a consistent approach. 
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Figure 7. Changes in separation with time between compo- 
nents C and J. Asterisks, diamonds, and triangles refer to 8.4- 
GHz VLBA data, 15-GHz VLBA data, and 4.8-GHz Space- VLBI 
(VSOP) data respectively. The solid line represents the regression 
fit to the 8.4-GHz VLBI data, while the dashed-lines represent the 
uncertainties from the fit parameters. Error bars are determined 
as in Section 3.3. 
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Figure 8. Changes in separation with time between com- 
ponents C and Nl, N2, and N3. The solid line repre- 
sents the regression fit to the data, while the dashed-lines 
represent the uncertainties from the fit parameters. Error 
bars are determined as in Section 3.3. Vertical lines show 



7-ray flares iPucella et al.l |2008|; ID'Ammando et all l2009al; 



CiDrini & Corbell I2009]: ID'Ammando et alj|2009bl: I 


-"ucella et al. 


2009: Vcrccllonc et al. 200sl ICutini & Havs 2009: 


Striani et al. 


2010). 



For each knot the time of zero separation obtained from 
the fit is T ,ni = 2008.61 ± 0.11, T , N2 = 2009.38 ± 0.17, 
and T ,N3 = 2009.93 ± 0.25 (Fig. |5J. Components Nl and 
N2 were also de t ected in VLBA observations at 43 GHz by 
iMarscher et al.l (|2010T ). Given the lower spatial resolution 
(and larger opacity) with respect to the 43 GHz, these 
components were detected at 15 GHz a few months later 
when their separation from the core was large enough to 
be resolved. The apparent separation speeds derived at 43 
GHz are 24±2 and 21.6±0.6 for Nl and N2 respectively, 
i.e. systematically larger than our values. This discrepancy 
may be due to the lower resolution of the 15-GHz data that 
m ay cause a b l ending with other features, as also suggested 
m lAbd o et al.l l|2010bf ). However, it must be noted that the 
higher separation speed found at 43 GHz is not reflected in 
a substantially different origin time of the jet components. 
In fact, in the case of Nl and N2 components the time of 
zero separation derived from 43-GH z data are 2008. 5 7i0.0 5 
and 2009.34T0.01 respectively (|Marscher et alj l201fj ). 
i.e. within the errors estimated from the 15-GHz data. 
In Table [3] we report the apparent separation speed and 
the ejection date of the superluminal knots of PKS 1510-089. 



Table 3. Proper motion of the jet components of PKS 1510-089. 
Column 1: source component; Col. 2: number of fitted epochs; 
Col. 3: apparent speed; Col. 4: zero-separation time; Col. 5: zero- 
scpara tion time in Julian D at e (JD); Col. 6: Refer enc e: 1: this pa- 
per: 2:lHoman et al] ll200lll: 3:lLister et al] d2009bl l: 4: lAbdo et all 
l|2010bh ; 5: IMarscher etail <201C|) . 



Comp 


N op 


/3app 


To 


To (JD) 


Ref. 


J 


8 


16.2±0.7 


1997.42i0.12 


2450602i45 


1 


A 


11 


15.0±0.7 


1994.92i0.18 


2449689i66 


1,2,3 


B 


19 


18.6±0.5 


2005.42i0.11 


2453524i40 


1,3 


Nl 


10 


17.3±0.9 


2008.61i0.11 


2454689i36 


1,4,5 


N2 


8 


18. Oil. 9 


2009.38i0.17 


2454971i62 


1,5 


N3 


5 


17.0i4.0 


2009.93i0.25 


2455171i91 


1 


C 


10 


20.2il.2 


1997.85i0.20 


2450759i73 


3 


D 


8 


18.9il.3 


1999.10i0.25 


2451215i91 


3 


E 


7 


14.8i0.9 


2000.20i0.20 


2451617i73 


3 


F 


6 


19.1il.5 


2004.42i0.28 


2453129i73 


3 



4 DISCUSSION 

The population of blazar objects represents a sub-class of 
AGN characterized by high levels of flux density variability 
across the entire electromagnetic spectrum. In the radio 
band the increase in luminosity is first detected at high fre- 
quencies, at millimeter wavelengths, subsequently followed 
at lower frequencies with some delay, consistent with opacity 
effects. Episodes of enhanced radio luminosity also seem to 
be related to changes in the pc-scale radio structure where 
new compon ents are ejected with apparent superl uminal 
velocity (e.g. IWagner et al.|[l995l : IJorstad et alj|200ll ). 



In the following we discuss our results on the multi-epoch 
analysis of the total intensity and polarimetric radio 
properties, and we compare the source activity at various 
frequency ranges. 



4.1 Superluminal motion and polarization analysis 

By means of the multi-epoch analysis of the parsec-scale 
structure of PKS 1510-089, various knots emerging from the 
core at different times have been detected. In particular, we 
obtain a solid estimate for one of the components detected 
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Figure 9. The 15-GHz lightcurve from VLBA MOJAVE data. 
The uncertainties on the flux density are within the symbols. Thin 
arrows repres ent the fit t ed ejec ted date of superluminal compo- 
nents from iLister et al.l ll2009bh , while thick arrows refer to the 
jet components studied in this paper. Each arrow refers to a dif- 
ferent superluminal component. The horizontal bars indicate the 
uncertainty on the ejection date. 



between 1999 and 2001, thanks to the addition of 5 new 
position measurements at 4.8 and 8.4 GHz. Interestingly, 
comparing the properties of the knots analysed in this 
paper with those studied in previous works (Table [3]) we 
found that all the knots are separating from the core with 
an apparent velocity that ranges between 15c and 20c, and 
all the jet components are well aligned along the same 
position angle, suggesting that no jet precession is taking 
place on a time lag longer than a decade in our frame. 
In Fig. [9] we report the time of zero separation of all the 
components ejected between 1994 and 2010, either analysed 
in this paper, or from the literature. The ejection of a new 
component likely occurs roughly once per year. The gaps 
between 2000 and 2004 is likely due to sparse observations 
available, precluding to reliably follow the evolution of the 
source components in this period. 

Since 1995, in addition to the typical variability interchang- 
ing low-activity with high-activity states, the source shows 
on average a slight increase of the total flux density (Fig. 
[9}. However, the sparse distribution of the data points does 
not allow us to unambiguously relate the ejection date of 
the superluminal components, represented by the arrows in 
Fig. [51 with the various states of the source variability. 
From the analysis of the parsec-scale resolution data 
we could construct the lightcurves for the core and jet 
components separately, and to study a connection between 
the total intensity and the polarized emission. The core has 
high degree of variability, both in the total intensity flux 
density and in polarization properties, but the data points 
are too sparse to compare changes in the total intensity 
and polarized emission throughout the entire time interval. 
An interesting case is represented by the simultaneous 
enhancement of the total intensity flux density and the 
decrease of the fractional polarization observed in the 
second half of 1999, and the rotation of EVPA around 
May- July 2001, with a swing of 85° and 130° at 8.4 and 
15 GHz respectively (Fig. [6]). Possible explanations may be 



related either to variation in the opacity in a newly ejected 
self-absorbed component, or to a highly ordered magnetic 
field produced by the compression of tangled magnetic 
fields by shocks. In the former scenario the evolution of the 
jet knot implies a change between the opacity regimes. As a 
consequence both the total intensity and the polarized flux 
densities should decrease during the transition, whereas the 
magnetic field should jump by 90 degrees. As the opacity 
decreases, the total flux density should increase in the 
optically-thick part of the spectrum, while the polarization 
percentage decreases. As the emitting region expands 
becoming optically-thin, the radio emission decreases, while 
the fractional polarization increases and the magnetic 
field flips of 90 degrees, going back to its original value. 
It must be noted that changes of the EVPA of 90° are 
expected from opacity effects during the transition between 
the two regimes. Since the process is quite fast, a dense 
time sampling of the polarization information (as the 
one obtained between 2007 and 2010, Fig. 1101) is crucial 
to detect the transition between the regimes. A sparser 
sampling, like the one between 1999 and 2001 (Fig. O, may 
provide EVPA changes larger/smaller than the expected 
90°. The transition between the optically-thick and -thin 
regimes occurs in the presence of high opacity values (e.g. 
r ~ 5 - 10, see for example Aller (1970) and Pacholczyk 
(1970) for a detailed discussion). Such opacity would cause 
a dramatic drop of the total intensity flux density during 
the transition. The lack of observations during the first 
half of 1999, i.e. when the new jet component should have 
originated, do not allow us to unambiguously confirm the 
opacity scenario. However, the increase of the total flux 
density observed between 1999 and 2000 is difficult to 
explain in the presence of such high values of the opacity, 
making this scenario unlikely. 

The alternative scenario assumes the formation of a 
shock that causes the compression of the plasma along 
the propagation axis. This generates the amplification of 
the perpendicular component of the magnetic field with 
respect to the parallel one, and an enhancement of the 
luminosity. This scenario would predict an increase of 
the fractional polarization instead of the drop observed 
close to the radio outburst. However, we must note that 
if the shock is oblique instead of transverse, the expected 
variations in the polarization properties are different, and 
strongly related to the obliqueness of the shock itself, and 
to the characteristics of the underling magnetic field like 
its order and strength. Therefore, a reliable interpretation 
of the observed polarization trends requires a much better 
sampling than that available between 1998-2001 (Fig. [6]), 
leaving the debate on the nature of the main mechanism at 
work still open. 

The analysis of the lightcurves and polarization trends 
between 1998 and 2001, shown in Fig. [Bl suggests that 
the variations of the total intensity flux density and the 
polarization properties may be explained by both the 
previous scenarios. Indeed, the total intensity and polar- 
ization properties can be related to the evolution (likely 
adiabatic expansion) of either a new jet component or a 
shock originated at the beginning of 1999. Support to this 
interpretation comes from the detection of a superlumi- 
nal component tha t possibly originated in 1999.10T0.25 
(|Lister et al l I2009bl 1. The rotation of the angle found in 
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2001 may be explained considering that this new component 
moves far enough from the core to be resolved and the 
intrinsic polarized emission of the core can be separated 
from that of the jet. 



4.2 Multifrequency analysis 

The lightcurves considered in the above explanation would 
require a more frequent sampling and for this reason a defini- 
tive interpretation of the physical mechanisms at work is 
precluded. A similar consideration may apply to explain the 
source behaviour aroun d April 2009, when information at 
other wavelengths (e.g. lAbdo et al.ll20 10b; Mar scher et al.l 
l20igk is availab l e. The m ultiwavelength lightcurve presented 
by lAbdo et all l|2010bh shows a flux density enhancement 
around that period that is first detected at 230 GHz, and 
after some delay at lower frequencies. The same behaviour 
is visible in Fig. [TT]where the lightcurves at 15 and 43 GH^f] 
have been compared. From Fig.[Tl]we see that the flux densi- 
ties at both frequencies have a similar behaviour, where the 
43-GHz data points seem to anticipate those at lower fre- 
quency. Interestingly, just after March 2009, when severa l 
7-ray flares have been detected llD'Ammando et alJ 1201 J ) . 
the flux density at 43 GHz becomes higher than that at 
15 GHz implying a change in the opacity of the compo- 
nent. The maximum of the radio emission is then reached 
in April 2009 when a nother strong 7-ray flare was detected 
l|Cutini fc H avs 2009). In correspondence of this luminosity 
enhancement the polarization percentage drops while the 
EVPA changes by about 75° (Fig. [TDJ at both 15 and 43 
GHz becoming parallel to the jet direction. In the same 
period a strong rotat ion of the optical polarization vector 
|Marscher et al.ll20ld) has been detected together with the 
ejection of a new superluminal comp onent observed at 1 5 
GHz (Fig. El), a nd already reported b v lAbdo et all (|2010bh . 
and at 43 GHz (jMarscher et alj|2010h . A similar example of 
dramatic rotation of the polarization angle and a drop of the 
fractional polarization in coincidence with a 7-ray flare was 
found in the blazar 3C 279, indica ting a co-spatiality of the 
optical and 7-ray emission region (jAbdo et alJl2010ch . This 
has been explained assuming a non-axisymmetric structure 
of the emitting region, likely a curved jet trajectory, rather 
than a perpendicular shock moving in an axially symmetric 
jet as we suggested in Section 4.1. In this case, the degree 
and angle of polarization strictly depends on the instan- 
taneous angle fo rmed by the direc tion of the motion with 
our line of sight (|Abdo et al.ll2010ch . The similarity between 
3C 279 and PKS 1510-089 suggests a common origin also for 
the emission at high energy and in the radio band, and the 
time lag in the flux density behaviour may be due to opacity 
effects as the shock passes by. A cospatiality of the 7-ray and 
radio emitting region was also claimed for the BL Lac ob- 
ject OJ 287, where two 7-ray flaring episodes o ccurred close 
in tim e with two major millimeter outbursts (|Agudo et al.l 
[2011]). 



3 All the values at 43 GHz presented in this paper are from 
the Large VLBA Project: Total & Polarized Intensity Images 
of Gamma-Ray Bright Blazars at 43 GHz, and available at 
[http://www.bu.edu/blazars/VLBA-GLAST/1510.htmll 
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Figure 10. Multi-epoch 15-GHz VLBA peak flux density, re- 
stored with the same beam size (upper panel), fractional po- 
larization (middle panel), and polarization angle (lower panel) 
of the core component of PKS 1510-089. Uncertainties are 
within the sy mbols. Vertical lines are relative to high-state 
7-ray activity jPucell a et al.l 120081: iD 'Ammando et al.l 1 2009a ; 
ICiprini fc C orbel 2009; D'Ammando et al.ll2009ri IPucella et al ' 
| 2009 l: IVercellone et al.l 120091 : ICutini fe Havd 120091 : IStriani et al 
I2OIOI) . 



In PKS 1510-089, the variation in the polarization angle at 
both 15 and 43 GHz suggests a change in the magnetic field 
orientation in the compact component rather than Faraday 
effects caused by an external screen. Changes in the po- 
larization angle of the sa me magnitud e at dif ferent wave- 
lengths were reported by iHoman et all l|2002bh who moni- 
tored the behaviour of a sample of 12 blazars by means of 
dual-frequency VLBA observations. 

It is worth noting that, as in the case of 3C 279 where no 
changes i n the radio band co uld be found related to the 7- 
ray flare jAbdo et all 12010(f ). in PKS 1510-089 no obvious 
conne ction between the 7-ray activity detected in March 
2009 ([D'Ammando et al.ll-01ll ) and the radio flux density 
behaviour has been found, suggesting that during this flare 
the synchrotron radiation in the radio band is not yet fully 
optically thin. 



5 CONCLUSIONS 

We have presented results from the analysis of multi-epoch 
polarimetric VLBI, Space- VLBI and archival VLBA data 
from the MOJAVE programme of the flat spectrum radio 
quasar PKS 1510-089 spanning over 15 years (1995-2010). 
This source shows a pc-scale core-jet structure where 
superluminal knots are ejected at different times. From the 
multi-epoch observations we found that the emission of 
new blobs is roughly constant with a time lag of about one 
year. Furthermore, the various jet components are moving 
away from the core with an apparent superluminal speed in 
the range between 15c and 20c and roughly with the same 
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Figure 11. Multi-epoch peak flux density at 15 GHz 
(squares) and 43 GHz (crosses) of the central component 
of PKS 1510-0 89. Vertical l i nes a re relative to high-state 7- 
ray activity llPucella et alj l2008t ID'Ammando et all l2009al: 
Ciprini fc Corbelll2009l: ID'Ammando et al.ll2009bl: IPucella et al. 
20091; IVercellone et aT]|2009l ; ICutini &: Hays! |2009|; IStriani et al. 



position angle, suggesting that the precession of the jet is 
not relevant on the timescale of decades in our frame. Both 
the total intensity and the polarized flux density of the 
core component show high level of variability. Our analysis 
shows that occasionally the EVPA has abrupt changes of 
about 90 degrees becoming roughly perpendicular to the 
jet direction during episodes related to enhanced radio 
emission and opacity variations. These properties may be 
explained assuming a change between the optically-thick 
and optically-thin regime as a consequence of a shock that 
varies the opacity. In this context, the luminosity locally 
increases due to the compression of the plasma in the di- 
rection of the shock propagation, causing the amplification 
of the perpendicular component of the magnetic field with 
respect to the parallel one, and the EVPA becomes parallel 
to the direction of the shock propagation. The observed 
properties may also be explained as due to a highly ordered 
magnetic field produced by either an oblique shock or 
a transverse one propagating down a jet with a curved 
trajectory. In this case the angle of the electric vector and 
level of polarization strictly depend on the instantaneous 
angle formed by the direction of the shock with our line 
of sight, thus explaining changes in the polarization angle 
different from the 90 degrees expected during the transition 
between the opacity regimes. However, the aforementioned 
scenarios cannot describe the source behaviour after all 
the detected flares and more monitoring multiwavelength 
campaigns are needed to properly understand the physical 
processes occurring in this source. 



APPENDIX A: MOJAVE MULTI-EPOCH 
MONITORING (1995-2010) 

To better characterize the evolution of the source structure 
with a unique approach we re-analysed VLBA data at 15 
GHz from t h e MOJ AVE programme al r eady published by 
iLister et ail l|2009rj ) and iHoman et all l|200ll) . Data from 
each epoch were modelfitted as described in Section 3.3, 
providing a homogeneous set of observations with the same 
data analysis procedure. 

The identification and monitoring of the source components 
is critical when two subsequent epochs are separated by a 
long time interval. Indeed, if the observation time coverage 
is sparse, it becomes very difficult to identify the same 
component at the various epochs. By comparing the multi- 
epoch visibility data we found that the source components 
could be reliably followed only in the observations from 
1995 July to 1998 December, and since 2007. 
Between 1995 and 1998 we could identify and follow a 
jet knot for which we derive an angular separation rate 
of 0.918±0.041 mas/yr, corresponding to an apparent 
separation ve l ocity /3 app = 15.0 ± 0.7. For this component 
IHoman eT al. (2001) derived an apparent speed of 15.7±0.4 
if converted in the cosmology used in this paper. From the 
linear back-extrapolation fit (Fig. IA1|) we estimate that the 
jet component was originated in 1994.92±0.18. 
Between 2007 and 2010, the pc-scale radio morphology dis- 
played a core-jet structure with the presence of four knots, 
three of which have been discussed in Section 3.3. The other 
knot we found increases its separation to the core with a 
rate of 1.134±0.035, corresponding to an apparent velocity 
/3app = 18.6 ± 0.5, and it should have been originated in 
2005.42 ±0.11 (Fig. ED). For this component ILister et afl 
( 2009b]) derived a separation speed of 17.0±2.6, and the 
ejection date 2005.87±0.25. The kinematic properties of the 
jet components are reported in Table [3] 
Between 1999 and 2006 the observations are too sparse in 
time to allow a reliable identification of the same source 
component throughout the various data sets. 
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